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ABSTRACT 

Using ray tracing for simple analytic profiles, we demonstrate that the lensing cross section for 
producing giant arcs has distinct contributions due to arcs formed through image distortion only, 
and arcs form from the merging of two or three images. We investigate the dependence of each of 
these contributions on halo ellipticity and on the slope of the density profile, and demonstrate that 
at fixed Einstein radius, the lensing cross section increases as the halo profile becomes steeper. We 
then compare simulations with and without baryonic cooling of the same cluster for a sample of six 
clusters, and demonstrate that cooling can increase the overall abundance of giant arcs by factors of a 
few. The net boost to the lensing probability for individual clusters is mass dependent, and can lower 
the effective low mass limit of lensing clusters. This last effect can potentially increase the number 
of lensing clusters by an extra 50%. While the magnitude of these effects may be overestimated due 
to the well known ovcrcooling problem in simulations, it is evident that baryonic cooling has a non- 
negligible impact on the expected abundance of giant arcs, and hence cosmological constraints from 
giant arc abundances may be subject to large systematic errors. 
Subject headings: cosmology: theory - lensing - giant arcs 
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1. INTRODUCTION 

Soon after the dis covery of giant arcs 
fl.vnris fc Petrosianl dH ISoucail et a,l.l \m& and 
their su bsequent interpr et ation as gravitational l ensing 
events (IPaczvnskilll987l) . |G rossman fc Naravanl <|198cf 
realized that the statistics of these relatively rare 
events can be used as a cos mological probe . This 
idea was greatly expanded by iMiralda-Escudd l)1993aj) 
and IWu fc Hammerl <) 19931) who used arc statistics to 



constrain the mass distributio n of clusters and the 
high redshift galaxy p opulation (|Miralda-Escudeiri993rJ: 
iBezecourt et al.lli998|) . Attempts to derive cosmological 
parameters from arc abundances soon followed, and 
concluded that giant arc abu ndances had a re latively 
mild cosmological dependence l)Wu fc Madl 199(f) . 

These early studies used simple analytical models of 
cluster mass distributions, which turn out to severely 
under-predict lensing prob abilities relative to c l usters 
in numerical simulatio ns ( B artelmann fc Weiss! 119941 
iBartelmann et aH T99iTL Nevertheless, it seemed reason- 
able to expect that the scaling of arc abundances with 
cosmology should be broadly consistent with predictions 
from analytic prescriptions, so it was surprising when the 
first numerical stud y to investigate arc abu ndances in 
various cosmologies l)Bartelmann et aTlll998j) found not 
only that there were order of magnitude differences in the 
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predictions from different cosmologies, but also that the 
abundance of giant arcs expected in the now-standard 
ACDM cosmology is an order of magnitude l ower than 
the observed abundance (|Le Fevre et al.lll994^ . 

The claimed sensitivity of giant arc abundances 
to cosmology and the strong inco nsistency between 
theory and observations found by IBartelmann et all 
( 1998) generated an explosion of both theoretical and 
observational studies. On the observational front, 
new arc samples with more cl early defined selec- 
tion functions became a vailable iLupnino et all 119991 
iZaritskv fc Gonzalezll2003l iGladders et alJl2003|K At the 
same time, the cosmological community mounted an ex- 
tensive theoretical effort to search for possible systemat- 
ics and/or theoretical uncertainties in order to determine 
whether the disagreement was real. 

The theoretical program has revealed that the problem 
of giant arc statistics is extremely rich, with predictions 
being sensitive to details of both background galaxies and 
cluster properties. Concerning sources, many models for 
giant arc abundances have assumed circular sources, but 
source ellipticity can s ignificantly boost the number of 
giant arcs in the sky l)Bartelmann et all 119951 iKeetonl 
|2001[) . Assumptions about th e distribution of source 
radii appea r to be impo rtant (|Miralda-Escudel Il993bt 
IQgurill200l IHo fc WhitdlSool . Th e redshift distribu- 
tion of sources plays a sig nificant role (jWambsganss et all 
120041 iDalal et alJ 120041) and in fact contributes much 
of the theoretical u ncertainty in modern predictions 
ijHennawi et al.ll2005|) . 

There are also several systematics associated with 
the lens population. For instance, while many studies 
have assumed spher i cally s ymmetric mass distributions 
e.g.. lHattori et all Il997t iHamana fc Futamasd 119971 



Molikawa et a,1J 11 99flt ICooravl 11 999alhl: IWilliams et al] 

1999), many others have pointed out that halo ellip- 



ticity significantly inc r eases gi ant arc abundance s (e.g. 
iGrossman fc Naravanl 119881: IBartelmann et alJ 119951: 
I l2002t iMeneghetti et aLll2003bl) . Moreover, halo 
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triaxiality can give rise to order of magnitude variations 
in the lensing probability of a given cl uster with viewing 
angle (jQguri et al J l2003t iDalal et al.ll2004f) . and there- 
fore must be treated properly in arc abundance predic- 
tions. Lensing probabilities for individual clusters ap- 
pear to be significantly enhanced during mergers (e.g., 
Torri et all l200i iMeneghetti etaD l2005bt iFedeli et all 



20061) . though such episodes appear to be rare and tran- 



sient enough that the overall lens populations is not sig- 
nificantly biased tow ards a cluster merging population 
ijHennawi et al.ll2005|) . 

One final related trend, noted in both analyti- 
cal and numerical works, is that lensing probabil- 
ities are very sensitive to t he radial density pro- 
files of lensi ng clusters (e.g., ^jrahia^sinjdej ll993a|; 
I Hattori et alJl!997t lOguri et alJl2001t IMeneghetti et all 
l2005aHHennawi et alJl2005|) . Consequently, clear predic- 
tions for the distribution of cluster density profiles must 
be in place for arc abundance predictions to become ro- 
bust. 

Systcmatics that do not seem to have significant effects 
on lensing probabilities have also been identified. For in- 
stance, studies in which a central galaxy is "painted" on 
the center of simulated clusters suggest that this addi- 
tional mass component has minimal impact on the ex- 



pected giant arc abundance l^ ojil^^^fc 
IMeneghetti et alJ2003aUDalal et al.l2004D . 
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presence of substructures and/or galaxies in the cluster 
also see ms to have a negligible effect on lensing prob- 
abilities llGrossman fc Naravanlll988t iFlores et al.ll2000t 
IMeneghetti et al.l l2000t | Hennawi et al. ||2005|). Line of 
sight projection effects llWambsganss et al.l2005D also ap- 
pear to be unimportant IjHennawi et al.ll2005l) . 
At present, there seems to be general agreement 

llDalal et all 12001 
120051). altho u gh th e 
situation is not completely clear <f Li et all 120051) . 
The new agreement arises from the coherent con- 
tributions of various effects, with no single system- 
atic expl aining the order of ma gnitude effects ob- 
served bv IBarte lmann et alJ lfl99ll). Recent analytica l 
ijKaufmann fc Straumann] l2000UBartelmann et al.ll2003|) 
and numerical <)Meneghetti et afll2005al) studies suggest 
that variations in the arc abundance between different 
cosmologies are factors of a few rather than an order 
of magnitude, further relaxing any real tension between 
theory and observations. What seems clear at this time 
is that arc abundances arc closer to old cosmology than 
to precision cosmology: agreement is to be understood 
at the factor of two level, even if current observations did 
not suffer from small number statistics. 

Given the complexity of the problem, it is important 
to consider whether there are any additional systematic 
effects that alter theoretical predictions. The goal of 
this paper is to study one such possibility: the impact 
of baryonic cooling on cluster mass profiles. Numerical 
simulations that include cooling have shown that once 
baryons have cooled and sunk to the center of the halo 
to form the cluster's central galaxy, the dark matter halo 
itself contracts in response to the increased central den- 
sity iBlumenthal et all 1198a iKochanek fc White! 120011 



T nerim et al . 112004 : ISellwood fc McCaughl \200K The 

effects of this contraction are significant, and include 
a steepening of the inner regions of the halo profile, a 



boost to the mass enclosed within a fixed radius r, and 
an increase of the degree of spherical sym metry of the 
mass distribution ( Ka zantzidis et all 12004^ . Since arc 
abundances are known to be sensitive to cluster den- 
sity profiles we expect that baryonic cooling may sig- 
nificantly affect the pr edictions. Indeed, recent work by 
iPuchwein et alJ l)2005[) has shown that baryonic cooling 
enhances the lensing probability of very massive clusters 
(M > 10 15 M ) by a factor of ~2-4. IPuchwein et"aT1 
l)2005[) did not consider, however, the effect of baryonic 
cooling on smaller clusters, which is particularly signifi- 
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In this work, we extend the analysis of the impact of 
baryonic cooling on arc statistics to clusters of consider- 
ably smaller masses. In particular, we demonstrate that 
baryonic cooling affects lensing probabilities for low mass 
clusters more than for high mass clusters, both because 
low mass clusters have a larger cooled mass fraction, and 
because lensing in low mass clusters is more confined to 
the central regions that are strongly affected by cooling. 
We also demonstrate that baryonic cooling has a strong 
isotropization effect for both the lensing probability of 
individual cluster as a function of viewing angle, and 
the spatial distribution of arcs around the cluster center. 
In order to understand why baryonic cooling can have 
such significant impact on the lensing cross section, we 
develop a theoretical framework that allows us to inter- 
pret the lensing probability associated with each cluster 
as a function of the minimum length to width ratio of 
the arcs considered. We have found that our theoretical 
framework can be used to better interpret and contextu- 
alizc results from previous works, providing an intuitive 
base upon which the sensitivity of arc statistics to var- 
ious cluster parameters can be deduced in a qualitative 
fashion. 

2. LENSING CROSS SECTIONS 

Lensing probabilities are usually characterized in terms 
of cross sections, so for completeness we review the con- 
cept of the lensing cross section. Let n arcs denote the 
number of giant arcs per unit area in the sky, expressed 
as, 

dzsdzid-p - — — — (1) 
dzsdzidp 

where dn aTCS / dz s dzidp is the number of giant arcs due 
to lenses at redshift zi, given sources at redshift z s with 
structure parameters p (such as radius and cllipticity). 
Our immediate goal is to compute c?n arcs /c?z s dz;dp. Let 
x denote a position vector on the source shell at red- 
shift z s behind a cluster at redshift z;. If r(x) is the 
length to width ratio of the image of a source at posi- 
tion x, then the expected number of arcs with a length 
to width ratio larger than £ behind this cluster can be 
written as dn arcs (£) = dn s a(£;p). Here dn s is the num- 
ber density of sources at redshift z s with parameters dp: 
dn a = (dn s /dp)(dV/dz s )dz s dp. Also, a(£\ p) is the area 
of the source plane over which the image of the source 
satisfies ?*(x) > £. This area o~(£) is called the lensing 
cross section, and is defined as 



a(£) 



r(x)>£ 



d 2 x ui(x; I) 



(2) 
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where w(x; £) is the number of images of a source placed 
at point x with a length to width ratio larger than £. 8 If 
we now denote the average lensing cross section of mass 
m halos as (a\m; z s , z\, p), it follows that the number of 
expected arcs is given by 



dn a 



dn s {p; z 8 ) dV 



dz s dzidp 



dp 



dz. 



(3) 



, dn-haios , \ dV i \ \ 

dm — {m,zi) — {a\m;z s ,zi 7 p) 7 

dm dzi 

where du^^ios/dm is the number density of halos of mass 
m at redshift z\. We can see that the lensing signal is 
fully characterized by the average lensing cross section 
(er|m; z s , zi,p) and the halo mass function dn/dm. 

To summarize, in order to compute the expected num- 
ber of giant arcs in the universe one needs to know: 

• The halo abundance dn/dm as a function of lens 
redshift z\. 

• The source abundance dn s /dp as a function of 
source redshift z s and structure parameters p. 

• The average lensing cross section for halos of mass 
m, as a function of lens and source redshift and 
source structure parameters. 

Given the evident complexity of the problem, studies 
that analyze the dependence of lensing cross sections on 
various parameters are particularly important: they help 
identify aspects of the theory that need to be well charac- 
terized to obtain robust predictions for arc abundances. 
As we shall see, the lensing cross section a(£) tends to 
exhibit some robust qualitative features, so in this pa- 
per we devote considerable effort to understanding these 
general features in order to provide a robust framework 
in which we can better interpret our numerical results. 

3. LENSING ALGORITHM AND SELECTION 
FUNCTION 

The idea behind our code is simple enough: simulation 
outputs arc used to construct three dimensional pixclizcd 
density maps of the clusters, which are then projected to 
create pixclizcd surface density maps. Using Fast Fourier 
Transforms, we compute the gravitational potential, an- 
gular deflection, and inverse magnificatio n tensor at each 
grid p oint in the lens plane. Following iHennawi et alJ 
l)2005[K we map each lens plane pixel back to the source 
plane, where we lay down a grid and associate each source 
plane pixel with all lens plane pixels contained within it. 
This yields a lookup table that allows us to quickly find 
which image (lens plane) pixels are lit given a set of lit 
source pixels. Circular sources are then placed along an 
additional coarser grid on the source plane, referred to 
as the placement grid. For each source, we identify its 
images and measure their lengths and widths. The lens- 
ing cross section for arcs above a given length to width 
ratio is then estimated from the number of placement 
pixels which generate appropriate arcs, weighted by the 
area of each placement grid pixel and the number of arcs 

8 Note that with this definition the phrase "lensing cross section" 
always refers to a cumulative quantity. The differential lensing 
cross section da/di is the multiplicity- weighted area over which 
sources generate arcs of length to width ratio exactly equal to I. 



produced by said source. All of the various parameters 
relevant for the algorithm (e.g., grid scales, grid extent, 
number of pixels, etc.) are chosen to ensure that the 
resulting cross sections are accurate to 10% or better. 
We note that because our goal is to estimate the rela- 
tive change in the lensing cross section of a cluster due 
to baryonic cooling, we do not expect our conclusions to 
be very sensitive to assumptions about source size and 
shape. We have therefore opted to make simple assump- 
tions, and leave a detailed analysis of how source prop- 
erties affect lensing cross sections to future work. 

The details of our algorithm are presented in Appendix 
lAl including various procedures to secure a 10% accuracy 
in the cross section, and several techniques to improve 
the speed of the algorithm. In this study, we have cho- 
sen to focus on tangential arcs only. This is not only 
because interpretation of the results is simpler when ne- 
glecting radial arcs (see Appendix [S] and the discussion 
in the following sections), but also because we expect the 
observational selection function for tangential arcs to be 
simpler than that of radial arcs since the former are more 
prominent and reside further away from a cluster's cen- 
tral galaxy than their radial counterparts. 

4. ANALYTIC MODELS 

We first study the cross section function a(£) for a va- 
riety of simple analytical models in order to develop a 
general understanding of its features. This basic frame- 
work will be used in section[5]to analyze the cross section 
function of numerical clusters. 

4.1. The Models 

We study halos with projected density profiles of the 
form 

_ b f(q)h((3) 



2 (q 2 x 2 + y 2 )^ 2 



(4) 



where q is the projected axis ratio, and b corresponds 
to the Einstein radius of the profile when q = 1. The 
function f(q) is defined by demanding that the mass con- 
tained within a circle of radius b be independent of the 
axis ratio q. The parameter (3 is the logarithmic slope 
of the profile, and h(f3) is defined by setting h(l) = 1 
and demanding that the mass contained within a circle 
of radius b be independent of (3. The special case (3 = 1 
corresponds to the Singular Isothermal Ellipsoid (SIE), 
and the spherically symmetric SIE case is called the Sin- 
gular Isothermal Sphere (SIS). For an SIS, the relation 
between Einstein radius and velocity dispersion is 



(5) 



where Di s and D os are the angular diameter distances 
between the lens and the source and between the observer 
and the source, respectively. 

In what follows, we consider two sets of models. 
For the first set, we fix (3 = 1 (SIE) and then con- 
sider Einstein radii corresponding to velocity dispersions 
ff„/(10 3 kms _1 ) G {0.6,0.7,0.8,0.9,1.0}. We also take 
q G {1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4} for a total of 35 models. 
For the second set, we fix the Einstein radius to that of an 
SIS with velocity dispersion a v = 10 3 km s _1 , and then 
consider density slopes (3 G {0.4,0.6,0.8,1.0,1.2,1.4}. 
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Fig. 1. — I map for a singular isothermal ellipsoid with cr„ = 
10 3 km s _1 and q = 0.8 aligned with the y axis. The map shows 
the function £(x), where £(x) is the largest I value of all tangential 
arcs created by a source centered on x. £ values are low for dark 
regions and high for bright regions. The white outside region was 
not searched for arcs. Finally, the dark, diamond shape solid line 
is the lens's caustic. There is a clear discontinuity in the i map 
all around the caustic, corresponding to sources which touch the 
caustic, and whose longest arcs arc therefore formed through the 
merging of two (fold arcs) or three (cusp arcs) distinct images. 
The dichotomy between distortion and image merging arcs is also 
reflected in the shape of the cross section curve cr(£), as seen Fig.l2l 

We use the same set of axis ratios, for a total of 35 ad- 
ditional models. In all cases we assume typical lens and 
source rcdshifts of zi = 0.3 and z s = 1.5. 

4.2. £ Maps and Cross- Sections 

Since the lensing cross section a{£) is simply the area 
of the source plane in which source produces arcs with 
length to widths larger than £, we begin our analysis by 
looking at £ maps for our analytic SIE models. 

The I map of a cluster is defined as the function £(x) 
on the source plane, where £(x) is the maximum length 
to width ratio among all tangential arcs produced by a 
source centered on x. Figure shows the £ map for an 
SIE profile with velocity dispersion o~ v = 10 3 km s" 1 and 
axis ratio q = 0.8. The topography of the map for the 
lens is quite striking. In particular, there is a discontinu- 
ity all around the lens's caustic (shown in the figure as a 
solid black line). The width of this "ribbon" is approxi- 
mately 2R where R is the radius of the source, indicating 
that all sources placed within this ribbon touch the caus- 
tic, and therefore correspond to sources in which distinct 
images merge to form a single long arc. 

We therefore introduce a distinction between arcs 
formed through distortion of a source, and arcs formed 
through the merging of several images. Moreover, a 
closer look at the ribbon around the lens's caustic reveals 
that there is an additional discontinuity around the lens's 
cusps, where not two but three images merge to produce 
a single arc. To differentiate between these two types of 
arcs, we shall refer to them as fold arcs and cusp arcs, 
respectively. 

The difference between distortion arcs and image merg- 
ing arcs is readily apparent in the shape of the cross 
section function a(£), shown in Figure [3 for various SIE 
profiles at fixed axis ratio q = 0.7 and for velocity disper- 
sions o-„/(10 3 km s _1 ) = 0.6,0.7,0.8,0.9, and 1.0. Each 
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Fig. 2. — Lensing cross section a(£) curves for SIE profiles with 
axis ratio q = 0.7. The velocity dispersions of the SIE halos are, 
from top to bottom, a v /(10 3 km s _1 ) = 1.0,0.9,0.8,0.7,0.6, and 
the assumed source radius is 5 kpc. The cross section function 
exhibits a transition between a roughly power law regime at low I 
values, and a roughly exponential regime at high I. The I value at 
which the transition occurs (It) corresponds to the scale at which 
the lensing cross section switches from being dominated by distor- 
tion arcs to being dominated by image merging arcs. The points 
show the actual results; the solid lines show phcnomcnological fits 
discussed in Appendix Ic^l 

cross section curve exhibits a clear transition, or knee, at 
some length to width ratio £t where it goes from being 
distortion dominated and roughly power-law, to image 
merging dominated and roughly exponential in nature. 
The value of this transition scale £ t clearly scales with 
the Einstein radius of the lens, reflecting the fact that 
sources will be more strongly distorted before touching 
the caustic for more massive clusters. 

Since the generic shape of the cross section function 
a{£) — a power law at low £, and a falling exponential at 
high £ — is mainly driven by the difference between dis- 
tortion arcs and image merging arcs, it is not surprising 
to learn that this shape is robust to changes in the halo 
parameters. Our main goal for the rest of this section 
is to understand how the various features of the lensing 
cross section function (its amplitude, the position of the 
transition scale, etc.) depend on the various halo param- 
eters, so that we may better interpret the cross section 
functions from the numerically simulated clusters. 

4.3. The Role of Ellipticity 
Figure shows the cross section function a(£) of three 
SIE profiles with velocity dispersion a v = 10 km s - 
and axis ratios q = 0.9 (diamonds), q = 0.7 (triangles), 
and q = 0.5 (squares). For illustration purposes, we have 
displaced the q = 0.9 data upwards by a factor of three, 
and the q = 0.5 data downwards by the same amount. 
The solid lines going through the points represent power 
law - exponential fits (see Appcndix[UJ). It is evident that 
the transition scale £< depends strongly on ellipticity. In 
particular, the more circular the profile, the larger the 
value for the transition scale £t- This reflects the fact 
that for a circular profile, merging images correspond 
to rings and typically have large £ ratios of order £ « 
(7r6)/(2i?). As the halo becomes more elliptical, arcs 
straighten out and the images no longer curve around 
the Einstein radius of the cluster, resulting in a smaller 
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Fig. 3. — Lensing cross sections for SIE profiles with a velocity 
dispersion a v = 10 3 km s _1 and axis ratios q = 0.5,0.7,0.9. For 
illustration purposes, we have displaced the q = 0.9 data upwards 
by a factor of three, and the q = 0.5 data downwards by the same 
amount. The solid lines are our best fit models. We see that the 
transition scale It moves towards lower I values as the ellipticity 
of the halo increases. Also shown with dotted lines are all cross 
section curves with their natural amplitudes. We see that as the 
halo becomes more elliptical, the amplitude of the image merging 
contribution to the cross section increases, reflecting the increase 
in the length of the tangential caustic. The amplitude of the image 
distortion contribution, however, appears to be roughly indepen- 
dent of halo ellipticity (at fixed Einstein radius). 

transition scale £±. 

Also shown in Figure [3] as dotted lines are the un- 
displaced (i.e., properly normalized) cross section curves 
for the q = 0.9 and q = 0.5 SIE halos. These curves 
demonstrate that as the halo becomes more elliptical the 
amplitude of the image merging contribution (high £) to 
the cross section increases reflecting the fact that more 
elliptical halos have longer tangential caustics. More sur- 
prising, however, is the fact that the amplitude of the im- 
age distortion contribution (low £) appears to be roughly 
independent of ellipticity, as evidenced by the fact that 
all three halos have almost identical cross sections at 
I w 5. We can understand this result qualitatively as 
follows: if we take a spherical halo and deform it into an 
elliptical one, the shear induced by the ellipticity will half 
the time add to the monopole shear, and half the time 
subtract from it. Consequently, the net cross section will 
be roughly unaffected, as illustrated above. 

There is one last point of interest that can be garnered 
from Figure [3] the maxim that elliptical profiles make 
better lenses is not always correct. More specifically, 
given a halo of a specified Einstein radius and a £ ratio 
of interest, the ellipticity with the highest cross section 
is that for which the transition scale £t is equal to the 
£ ratio of interest. For instance, since the cross section 
quickly falls beyond £ t , at high £ values (£ > nb/2R) and 
fixed Einstein radius only circular profiles are effective 
lenses. 

4.4. The Role of the Density Slope 

To understand how the slope of the density profile af- 
fects lensing, we first look at the £ maps of a series of 
analytic models. Figure 0] shows the maps for various 
power-law profiles where we fixed the Einstein radius to 
that of an SIS with velocity dispersion a v = 10 3 km s _1 




Fig. 4. — £-maps for power-law profiles with axis ratio q = 0.6, 
for profile slopes = f.4 (top-left), = 1.2 (top-right), = 1.0 
(middle-left), = 0.8 (middle-right), = 0.6 (bottom-left), and 
= 0.4 (bottom-right). The long axis of the lens is always along 
the y axis. The Einstein radius is fixed to that of an SIS with a 
velocity dispersion cr v = 10 3 km s . We have used the same gray 
scale in all maps to allow easy comparison between them. Note that 
fold arcs and distortion arcs not aligned with the long axis of the 
lens quickly become smaller (have smaller I values) as decreases. 
For extremely shallow profiles, this reduction in I can be so severe 
that fold arcs may become much smaller than single image arcs 
along the long axis of the lens, and the fold arc contribution to the 
cross section becomes unresolved. As a final note, the dark rings 
seen in some of the maps are due to sources placed along the radial 
caustic after we filter out the I values from the corresponding radial 
arcs. 

and the axis ratio to q = 0.6. In order to facilitate com- 
parison between the various maps, we have used the same 
gray scale in all panels. 

We see that the division of arcs into single image, fold, 
and cusp arcs is generic. What is striking is the strong 
impact of the profile slope on the length to width ratios 
of fold, and to a lesser extent, cusp arcs. The "ribbon" 
around the lens's caustic decreases monotonically in in- 
tensity as (3 decreases, and is even truncated in the case 
of the (3 = 0.4 profile. The cusp arcs of the flatter profiles 
also appear to have smaller length to width ratios than 
their counterparts in steep profiles, but the difference is 
less extreme. 

A similar pattern is also evident in the image distortion 
component of the maps, where we see that the length to 
width ratio of sources not aligned with the long axis of 
the lens decreases quickly as the profile becomes shal- 
lower. We can explain these patterns through a heuristic 
argument as follows. Consider a small source located off 
the long axis of the lens. If the halo is highly elliptical, 
the local curvature of the mass distribution is low and 
hence the shear induced through ellipticity is low. Like- 
wise, flattening of the profile reduces the shear induced 
by the mass distribution by making it more uniform, so 
flat, elliptical profiles are poor lenses if the source is off 
the long axis of the lens. If a source is along the long axis 
of the lens, however, the curvature of the mass distribu- 
tion is high and the resulting images are strongly dis- 
torted, even if the profile itself is shallow. Thus, highly 
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Fig. 5. — Cross sections for power law profiles of various slopes 
/3. The axis ratio is fixed at q = 0.6, and the Einstein radius is 
that of an SIS with a velocity dispersion a v = 10 3 km s —1 . The 
corresponding £ maps can be seen in Fig.|I] Note that the curves 
for f) > 0.8 are all similar in shape, with the transition scale it de- 
creasing monotonically with /3, reflecting the lower I value of fold 
arcs for shallow profiles. At /3 = 0.6 the profile becomes so flat that 
single image arcs along the long axis of the lens can be substantially 
larger than fold arcs. Consequently, the single image contribution 
to the cross section "spills over" into the fold arc regime, qualita- 
tively changing the features of the cross section. Further flattening 
of the profile — the f3 = 0.4 case — results in such small fold arcs 
that their contribution to the cross section becomes completely 
washed out by the single image contribution. The cusp arc contri- 
bution, however, becomes prominent, and the generic shape of the 
cross section is once again similar to that of the steeper /3 > 0.8 
profiles. 

elliptical flat profiles will only produce arcs when the 
sources are placed along the long axis of the lens. 

Let us now turn our attention to the cross section func- 
tion o~(£), shown in Figure [3] for the sample cases from 
Figure^] Focus first on the top four curves, correspond- 
ing to (3 = 1.4,1.2,1.0, and 0.8. All four of these have 
the characteristic two-component shape, with the transi- 
tion scale it decreasing as the profile becomes shallower. 
This agrees with our expectations: the transition scale It 
for these curves is roughly the minimum length to width 
ratio of fold arcs. Since fold arcs become smaller as the 
profile becomes shallower, the corresponding transition 
scale decreases. 

Turning to the (3 = 0.6 curve, however, we see that 
there is no obvious transition. Furthermore, the (3 = 0.4 
curve has a transition that occurs at larger values than 
the transition for the (3 = 0.8 case. This behavior may 
seem puzzling at first, but is easily understood as fol- 
lows. As j3 decreases from (3 = 0.8, fold arcs keep get- 
ting smaller and smaller, so at some point distortion arcs 
from sources placed along the long axis of the lens be- 
come longer than off-axis fold arcs. In other words, the 
distortion component of the cross section will "spill over" 
into the fold arc region, thereby erasing the clear tran- 
sition observed in the steeper profiles. As j3 is further 
decreased, fold arcs become negligible, and a new transi- 
tion becomes evident, but now the transition is between 
distortion arcs and cusp arcs. This explains why the 
= 0.4 transition occurs at higher length to width ra- 
tios than the (3 = 0.8 transition. 

There is an additional trend apparent in the cross sec- 
tion functions shown in Figure|S| as the profile gets shal- 



lower, the amplitude of the image distortion contribu- 
tion to the cross section decreases with f3, even though 
the Einstein radius of the various profiles has been held 
fixed. This reflects the fact that flat profiles induce lit- 
tle gravitational shear, and are thus not very effective at 
distorting the images of sources. 

5. SIMULATED CLUSTERS 
5.1. The Cosmological Simulations 

Next, we analyze high-resolution cosmological simula- 
tions of six cluster-size systems in the "concordance" flat 
ACDM model: fi M = 1 - A = 0.3, Q b = 0.04286, 
ft = 0.7 and <j$ = 0.9, where the Hubble constant is 
defined as Hq = 100 ft km s" 1 Mpc -1 , and erg is the 
power spectrum normalization on an 8ft -1 Mpc scale. 
The simulations were done with the Adap tive Refinement 
Tree (ART) A-body +gasdynamics code l)Kravtsovll999l 
iKravtsov et al.l20"f)2l) . a Eulerian code that uses adaptive 
refinement in s pace and time, and (non-adaptive) refine- 
ment in mass l|Klvpin et al.ll2001jl to reach the high dy- 
namic range required to resolve cores of halos formed in 
self-consistent cosmological simulations. The simulations 
presented here a re a subset of th e simu lated cluster sam- 
ple pr esented in IKravtsov et all l)2006|) and lNagai et alJ 
(2006), and we refer the reader to these papers for more 
details. Here we summarize the main parameters of the 
simulations and list the basic properties of clusters at 
z=0 in Table □ 

High- resolution simulations were run using a 128 3 uni- 
form grid and 8 levels of mesh refinement in the computa- 
tional boxes of 120ft," 1 Mpc for CL1-CL3 and 80ft," 1 Mpc 
for CL4-6, which corresponds to the dynamic range 
of 128 x 2 8 = 32768 and peak formal resolution of 
80/32768 w 2.44ft -1 kpc, corresponding to the actual 
resolution of sw 2 x 2.44 « 5ft -1 kpc. Only the region 
of ~3-10ft -1 Mpc around the cluster was adaptively re- 
fined, while the rest of the volume was followed on a uni- 
form 128 3 grid. The mass resolution corresponds to the 
effective 512 3 particles in the entire box, or the Nyquist 
wavelength of \-^ y = 0.469ft -1 and 0.312ft -1 comoving 
Mpc for CL1-3 and CL4-6, respectively, or 0.018ft -1 and 
0.006ft -1 Mpc in the physical units at the initial redshift 
of the simulations. The dark matter particle mass in the 
region around the cluster was 9.1 x 10 8 ft -1 M Q for CL1-3 
and 2.7 x 10 8 ft _1 M Q for CL4-6, while other regions were 
simulated with lower mass resolution. 

We repeated each cluster simulation with and without 
radiative cooling and processes of galaxy formation. The 
first set of "adiabatic" simulations have included only 
the standard gas dynamics for the baryonic component 
without gas cooling and star formation. The second set 
of simulations included gas dynamics and several phys- 
ical processes critical to various aspects of galaxy for- 
mation: star formation, metal enrichment and thermal 
feedback due to supernovae (type II and type la), self- 
consistent advection of metals, metallicity dependent ra- 
diative cooling an d UV heating due to cos mological ion- 
izing background (jHaardt fc Madaull 19961) . Throughout 
this paper, we refer to the adiabatic simulations and sim- 
ulations with cooling and star formation simply as "adi- 
abatic" and "cooling" run, respectively. 

In the cooling run, the cooling and heating rates take 
into account Compton heating and cooling of plasma, UV 
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TABLE 1 

Properties of simulated clusters from cooling runs 
at z=0 





-^500c 

(h- 1 Mpc) 


JW 500c 
(lO 13 /^ 1 Mq) 


» rtot 
JW 500c 

(lO 14 ^- 1 Mq) 


CL1 


1.160 


8.19 


9.08 


CL2 


0.976 


5.17 


5.39 


CL3 


0.711 


1.92 


2.09 


CL4 


0.609 


1.06 


1.31 


CL5 


0.661 


1.38 


1.68 


CL6 


0.624 


1.22 


1.41 



heating, and atomic and molecular cooling and are tab- 
ulated for the temperature range 10 2 < T < 10 9 K and a 
grid of metallicitics, and UV intensities u sing the Cloudy- 
code (ver. 96b4, iFerland et all fl998^ . The Cloudy 
cooling and heating rates take into account metallic- 
ity of the gas, which is calculated self-consistently in 
the simulation, so that the local cooling rates depend 
on the local metallicity of the gas. Star formation in 
these simulations was done using the observationally- 
motivated recipe fe.g.. lKennicuttJll998|) : p* = Pg a s/i*, 

with t» = 4 x 10 9 yrs. The code also accounts for the 
stellar feedback on the surrounding gas, including injec- 
tion of energy and heavy elements (metals) via stellar 
winds and supernovae and secular mass loss. 

Starting from the well-defined cosmological initial con- 
ditions, these simulations follow the formation of galaxy 
clusters and capture the dynamics of dark matter, stars 
and ICM self-consistently. These simulations can there- 
fore be used to examine the effects of ga s cooling and 
star fo rmation on the density profiles (e.g., IGnedin et alJ 
12004(1 and shapes l|Kazantzidis et al.ll2004|) of dark mat- 
ter halos as well as the roles of substructure present in a 
realistic cosmological context. Although the magnitude 
of the effects are likely overestimated in the current sim- 
ulations due to overcooling problem, the effects of baryon 
cooling we discuss below are generic. 



5.2. Surface Density Maps and Ray Tracing 

To create the surface density map of a cluster, we be- 
gin by nesting 2/i _1 Mpc and 8/i _1 Mpc boxes on the 
center of the cluster (defined as the position of the most 
bound dark matter particle in the halo). The boxes have 
256 and 64 pixels to a side, respectively, and the density 
at each grid point is computed with a cloud in cell al- 
gorithm. The 3D density maps are projected along the 
region < r/(l/i -1 Mpc) < 1 and 1 < r/(l/i _1 Mpc) < 4 
for the small and large boxes respectively, and the result- 
ing surface density maps added to produce a 2D surface 
density map which is 8/i _1 Mpc in extent and reaches a 
resolution of 7.8/i _1 kpc in the central regions. 

Given the surface density map, we use Fast Fourier 
Transforms to compute the angular deflection and in- 
verse magnification tensor at each lens plane grid-point. 
We then select a small region near the center of the 



lens plane, and refine the lens plane grid using bilin- 
ear interpolation. The refined lens plane resolution r\ 
is defined in terms of the source plane resolution r s via 
ri = ?* s /max(AR,) where r s is the source plane grid reso- 
lution and max(AR,) is the maximum value of the radial 
eigenvalue Ar over all pixels where |Ar/At| > 0.9£ m i„ 
(i.e., all pixels relevant for lensing). The source plane 
resolution r s is itself defined via r s = R/N where R is 
the source radius and N is the number of pixels that 
"fit" within R, for which we take N = 5^2 as our de- 
fault value. Finally, once the lens and source plane grids 
arc defined, we link each source plane pixel to every lens 
plane pixel that maps to it, resulting in a lookup table 
that can quickly generate the lensed image of any source 
on the source grid. 

5.3. The Impact of Gas Cooling in Lensing Cross 
Sections 

We start by analyzing one of the relaxed clusters with 
an isolated central density peak to illustrate the general 
features. Figure shows the £ maps of the adiabatic and 
cooling runs of CL2, viewed along a given line of sight. 
The difference between the two maps is dramatic: the 
adiabatic cluster is strongly asymmetrical, and arcs can 
only be produced if sources are placed along the long axis 
of the lens. The cooling cluster, by contrast, can produce 
arcs for sources at all azimuthal angles. The cross section 
for fold arcs is very small in the adiabatic cluster, but 
prominent in the cooling cluster. These are exactly the 
trends we would expect based on our analytic models 
given that adiabatic clusters are flat and highly elliptical, 
whereas cooling clusters tend to be less elliptical and have 
much steeper inner density profiles. 

Next, we study the effect of gas cooling and star forma- 
tion on the cross section function a(£). Figure shows 
the average cross section for both the adiabatic and cool- 
ing simulations. The impact of baryonic cooling is imme- 
diately evident: at I = 5, the boost in the cross section 
due to cooling is a factor of f=a 4. The cooling cluster is 
capable of forming much larger arcs than the adiabatic 
cluster because gas cooling enhances the mass in the cen- 
tral regions of the cluster, steepens its density profile, 
and makes the halo less elliptical. All of these effects 
make the transition scale £t shift to larger values. The 
net effect is a dramatic increase in the length to width 
ratio of the arcs that the cluster can produce. Turning 
this around, we infer that the minimum mass cluster that 
can produce arcs above some given £ is smaller with cool- 
ing than without. Reducing the effective mass cut may 
significantly increase the number of lensing clusters, de- 
pending on the steepness of the halo mass function. We 
make an estimate of this effect in section HJ 

Also striking in Figure [3 is the difference in the vari- 
ance in the cross section among different lines of sight: 
the cooling simulation exhibits considerably smaller vari- 
ance than the adiabatic one. In the absence of gas cool- 
ing, the lensing cr oss section varies strongly with view- 
ing angle (see also lDalal et al.lf2004tl . This is due to the 
combined effects of the flatness of the NFW profile and 

9 Operationally, a small square section of the lens plane is se- 
lected by demanding that no pixel within 50 kpc of the edge of the 
square have an eigenvalue ratio |Ajj/Ax| > 0.9£ m i„. This ensures 
all pixels relevant for lensing are well within the selected region of 
space. 
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Fig. 6. — £ maps for CL2 in the adiabatic (top) and cooling (bot- 
tom) simulations. We have used the same gray scale for both maps 
to facilitate comparison. The difference between the two maps is 
dramatic. In particular, adiabatic clusters can produce arcs only 
for sources located along the major axis of the lens, and these are 
cusp arcs; fold arcs provide a negligible contribution to the lensing 
cross section. Cooling clusters, on the other hand, exhibit a much 
more isotropic arc distribution, and should exhibit prominent fold 
arcs. These qualitative features are exactly what we would expect 
given that adiabatic clusters tend to be flat and highly elliptical, 
whereas cooling clusters are less elliptical and have much steeper 
profiles in the inner regions of the cluster. These maps are char- 
acteristic of all adiabatic and cooling maps we examined in which 
the central density peak was relatively isolated. 



triaxiality of a halo. That is, if the inner density profile 
is shallow, a small boost to the density due to changes in 
the viewing angle can dramatically increase the Einstein 
radius and hence the lensing cross section, resulting in a 
large dispersion in the cross sections for different lines of 
sight. In the cooling simulations, by contrast, the Ein- 
stein radius and cross section of the lens remain relatively 
unchanged as we vary the viewing angle. The main rea- 
son is that gas cooling steepens the inner density profile 
of the cluster, making the cross section less sensitive to 
boosts in the projected mass. 10 The fact that the cooling 
clusters are less triaxial also contributes to the reduced 
variance, but is not a dominant effect. 

10 Recall that at fixed slope and Einstein radius, the amplitude of 
the distortion componen t of t he lensing cross section is independent 
of cllipticity (see section [4.ijl . 




Fig. 7. — The average cross section and its variance among 
25 random lines of sight for a massive cluster with a relatively 
isolated central density peak in cooling (solid line) and adiabatic 
(dashed line) simulations. Cooling boosts the average cross section 
by a factor of 4 at low £ and significantly more at higher £. In 
addition, it significantly reduces the variance among lines of sight. 
Finally, note that the cooling cluster is a much more effective lens, 
capable of producing significantly larger arcs than the adiabatic 
cluster. This implies that cooling decreases the lower mass cut in 
lensing-selected clusters. 

5.4. The Impact of Substructures 

We now examine the effects of substructure on the lens- 
ing cross section. Figure |H] shows two typical £ maps 
for clusters with substructure, for both adiabatic (top 
panel) and cooling (bottom panel) simulations. In both 
cases, substructure leads to the appearance of "bridges" 
between the central density peak and the most nearby 
substructure. 

It is apparent from the maps, however, that such 
bridges provide only a modest contribution to the lensing 
cross section of the largest arcs in both the adiabatic and 
the cooling clusters, in qualitative agreement with previ- 
ous estimates of the impact of substructure on the lensing 
cross section (iMeneghetti et alJl2000l iFlores et al.ll2000i 
iHennawi et al.ll2005|K ~ 

Figurc[3]shows the average cross section for cluster CL1 
over 25 random line of sight projections. Even though 
the central density peak of CL1 is not isolated, the im- 
pact of baryonic cooling on the lensing cross section is 
strikingly similar to that of CL2, which did have an iso- 
lated central density peaks (cf. Fig.[7|). In particular, the 
lensing cross section is still boosted by a factor of w 3 
in the cooling run, and the transition length to width 
ratio £ t is again shifted to larger values. We found this 
to be the case in all of our simulated clusters. We con- 
clude that substructures do not significantly change how 
baryonic cooling affects arc formation. A complete quan- 
titative characterization of these effect, however, requires 
a larger cluster sample than what is available to us at this 
time. 

6. THE IMPACT OF BARYONIC COOLING ON 
GIANT ARC ABUNDANCE ESTIMATES 

As discussed in scction|21 estimating the number of gi- 
ant arcs per unit area in the sky involves characterizing 
the average lensing cross section (a\m; zi, z s , p) for halos 
of mass m as a function of rcdshift and for a range of 



Impact of Baryonic Cooling on Giant Arcs 



9 



Adiabatic Simulation 



10" 



10" 



10" 



Cooling run 
Adiabatic run 




-100 o 



Cooling Simulatio 




Fig. 8. — I maps of clusters in which there is substructure near 
the main density peak. The top panel shows an adiabatic clus- 
ter, and the bottom panel shows a cooling cluster (not the same 
cluster). The basic map topology around the central density peak 
is largely unaffected by the substructure, though we often found 
a "bridge" connecting the central density peak with the nearest 
substructure. This bridge enhances the lensing cross section for 
moderate length to width ratios, but in general we expect the ef- 
fect to be modest. 

source structure parameters (e.g., radius and ellipticity). 
Since our lens sample is limited to a handful of clusters, 
we cannot provide a robust estimate of (a\m; zi, z s , p), 
even ignoring dependencies on source structure parame- 
ters. It is thus evident that extrapolation of our results 
is necessary, and will involve some guesswork. In light 
of these difficulties, we have opted to use the simplest 
possible scaling arguments to estimate the ratio of the 
number of giant arcs expected in cooling simulations to 
the number expected in simulations with only adiabatic 
gas physics. 

Our analysis begins by noting that eq. © can be writ- 
ten as a product, 

dn nrrx 

g{z s ,zi,p) T{£;z s ,z h p), 



where 



and 



dz s dzidp 



g(z s ,z h p) 



(6) 



dn s (p;z s ) dV dV 



dp 



dzc. dz 



(7) 



■ ZhP) 



dm — — <r(i\m, z s 
dm 



zuP)- (8) 



Fig. 9. — The average cross section and its variance for our most 
massive cluster (CL1). This cluster has substructures present near 
its central density peak, yet the average cross section function a(£) 
exhibits the same general features seen for the isolated central den- 
sity peak case examined in Fig. |7| We found this to be true of all 
the clusters we considered. Thus, the presence of substructures 
does not appear to significantly affect the impact of baryonic cool- 
ing on arc abundances. 



The quantity r is the lensing optical depth, and it con- 
tains all of the effects of baryonic cooling on arc abun- 
dances. (The function g(z s , Zi, p) does not depend on the 
internal structure of the lenses.) In particular, the ratio 
of arcs in the cooling and adiabatic simulations is simply 



N ad 



Tad 



oi _ J dm (dn/dm) a coo i(l\m) 
J dm (dn/dm) a a d((\m) 



(9) 



Our goal is to estimate this ratio. Even this simple sce- 
nario requires significant extrapolation from our small 
sample of clusters, so we again opt for making the argu- 
ment as simple as possible. We have then 



Tad 



dm -P-a a d(£\m) 
dm 



(10) 



where m a d is the mass cutoff for the adiabatic case and 
we assume a power law scaling a a d — C a dmP '■ For an 
SIS lens, 7 ~ 1.3. For a power law halo mass function 
dn/dm = Am~ a we have then 



Tad 



ACad 



1 



a — 7 - 1 rr 



a— 7 — 1 ' 
ad 



(11) 



A similar expression holds for t coo ;, with C coo i = XC a d 
where A is the net boost to the lensing cross section due 
to cooling, A w 3. Wc see then that 



Tad 



A 



m ad 

Wlcool 



a— 7 — 1 



(12) 



The second factor in parenthesis on the right hand side 
represents the additional boost to the optical depth due 
to the lowering of the halo mass cutoff. To give a rough 
estimate for this term, note that the cutoff mass must 
occur when I ks i t (m). We know that £ c t ool (m) > if(m), 
but due to our small cluster sample, we do not know how 
it scales with mass in either case. However, we can make 
some progress if we assume a power law scaling between 
it and m, and further assume that cooling only changes 
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the amplitude of this scaling by some fixed factor f. 11 
In that case, the mass cutoffs m a d and m coo i are related 
via 

km P ad = i = fkm cool (13) 

where f3 is the slope of the relation £t(m) = km" for 
adiabatic clusters. For spherical halos, we know that 
£ t w nb/2R, so in general we expect l t oc b cx a 1 / 2 oc m 7 , 
and hence (3 w 7/2. Finally, the ratio / = l c t ool jt a t d ss 1.5 
in the simulations, which results in 

Tcooi _ ^2(a-7-X)/7_ 
Tad 

For a k -3 as appropriate for clusters, and 7 ~ 1.3 as 
expected for SIS lenses 12 , we find (a — 7 — 1)//? w 1, 
meaning the decrease in the halo mass cut due to cool- 
ing can further enhance the number of arcs in the sky 
by about 50%. This is smaller than the net increase due 
to the boost in cross section, but still far from negligi- 
ble. Moreover, for extremely large arcs, the correspond- 
ing mass scale will be larger and hence the mass function 
will be steeper, thereby enhancing the importance of the 
change in the mass cutoff scale due to baryonic cooling. 

7. DISCUSSION AND COMPARISON WITH 
PREVIOUS WORK 

Perhaps the most obvi ous point of comparis on for 
our work is the work of iPuchwein et al.1 l|2005|) , who 
studied the impact of baryonic cooling on very massive 
(M > 10 15 M Q ) clusters. To the extent that our analy- 
sis overlap, our results are in agreement and we concur 
with their conclusions: baryonic cooling boosts the lens- 
ing cross section of the most massive clusters by a factor 
of ?s 2 — 4, and this boost is primarily driven by the steep- 
ening of the profile. For smaller mass clusters, this boost 
can be significantly larger, which can in turn significantly 
increase the number of arcs systems in the universe de- 
pending on the minimum length to width ratio of the 
arcs considered. 

In addition to the work by IPuchwein et alJ l)2005fl . 
there have been several other studies that investigated 
the impact of the presence of the mass ive central galaxy 
found at the center of most clusters l |M eneghetti et alJ 
l2003bllDa,lal et al J2004tlHo"fc Wbitell200^. These stud- 
ies noted that when a mass model for the central cluster 
galaxies is tacked on numerically simulated dark matter 
halos the lensing cross section for the cluster is relatively 
unchanged. This demonstrates that the halo's response 
to the baryonic concentration at the center is critically 
important, and strongly supports our argument that the 
enhancement to the lensing cross section is driven by the 
associated steepening of the halo profile. 

Incidentally, we expect that the steepening of the 
halo profile induced by cooling will have an additional 
interesting effect which we have thus far ignored: it 
should increase the sensitivity of lensing cross sect i ons to 
the assumed source redshift. iWambsganss et all (|2004fl 
used the approximation [i w |Ar/At| ~ i in esti- 
mating the impact of source redshift on lensing cross 

11 Note that since we do not know the scaling of it with mass in 
either case, this assumption represents the simplest possible rela- 
tion between the two scalings, and is not guaranteed to be correct. 

12 a oc b 2 oc at oc m 4 / 3 for an SIS lens 



sections, and demonstrated that a broad redshift dis- 
tribution enhan ces the lensing optical depth. Subse- 
quent work Ce.g.. lDalal et al.l2004tEi et all2005|) demon- 
strated th at the sensitivity to source redshift was overes- 
timated in lWambsganss et al.l ((2004ft since their approxi- 
mation fi ?s £ is exact only for SIS profiles, with shallower 
profiles resulting in less sensitivity to the source redshift 
distribution. Given that baryonic cooling tends to both 
steepen and circularize the halo mass profile, we expect 
the cooling clusters to be more sensitive to the assumed 
redshift distribution. 

Another interesting way in which our work can be 
related to previous studies is to ask whether the gen- 
eral theoretical framework developed here can shed some 
light on previous results. For instance, we have already 
seen that the question of the isotropy of the lensing 
cross section is a strong function of the sharpness of the 
ha lo density profile, in agreement with the arguments 
bv lDalal et alJ l)2004t) for flat, dark matter only profiles. 
Likewise, our theoretical framework can help explain how 
the lensing cross section scales with mass. At low masses, 
the arc cross section is dominated by merging arcs, and 
thus should quickly drop with decreasing mass. Con- 
versely, above the mass cutoff, the lensing cross section 
should become dominated by distortion arcs and have a 
much more mild mass dependence. This characteristic 
behavior of a rapidly rising cross section at low masses 
transitioning into a flatter regime seems to be at least in 
qualitati ve agreement wi t h the cross section mass scaling 
found bv lHennawi et all i|2005ft . 

This same type of reasoning can be used to estimate 
the dependence of lensing cross section on source radii. 
For instance, consider first the case in which the mini- 
mum length to width ratio of interest £ <^ £ t where £t is 
the transition scale for a given cluster and source size. In 
this limit, the cross section is completely distortion dom- 
inated and hence source size independent (sec appendix 
IbI for details) . If I < it , increasing £ t will slightly boost 
the image merging cross section, so we expect the lensing 
cross section to increase with increasing source size. Con- 
versely, at i > it, the length to width ratio of a source 
at a fixed distance d < R from the caustic will decrease 
with increasing source radius due since the source then 
extending out to regions of lower eigenvalue ratio. This 
implies that the contour of fixed £ moves closer to the 
caustic, and hence the lensing cross section will decrease 
with increasing source radius. 

Based on these arguments, we expect the scaling of 
the lensing cross section <r{£) with source radius to 
be quite complicated. This brings up the question 
of why should we only consider arcs as a function of 
their length to width ratio? While early works (e.g., 
iGrossman &; Naravanlll98cllMiralda-Escudjll993aj) con- 
sidered arc abundances as a function of arc length, width, 
and length to width ratio, this is no longer the case. We 
expect that the move towards considering only length to 
width ratios came from the naive expectation that the 
correspondin g cross section would b e source radius inde- 
pendent fe.g. lWu fc Hammerlll993|) . Such a result would 
imply that arc abundances could be predicted without 
detailed knowledge of the radius distribution of high red- 
shift galaxies, an obscrvationally challenging problem. 
Unfortunately, <j{£) does depend on source radii, so the 
advantage of considering a(£) only is considerably less- 
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Fig. 10. — The lensing cross section for an SIE profile with axis 
ratio q = 0.8 is shown above with diamonds. Also shown as trian- 
gles is our estimate for the distortion contribution to the lensing 
cross section obtained from our best fit model to the data (see Ap- 
pendix |C}. The solid line is ob tained using the semi-analytic pre- 
scription of lFcdcli ct al.l 120061) . while the dashed line is the source 
plane area over which the eigenvalue ratio satisfies |Aj^/Ax| > I- 
Finally, the dotted line is the result from the semi-analytic pre- 
scription of LFcdel i ct al.l 120061) after subtracting the amplitude at 
the transition scale It as determined by our best fit. As expected, 
the resulting cross section is in excellent agreement with our best fit 
estimate of the distortion contribution to the lensing cross section. 



ened. Rather, given that the source radius dependence 
of o~(L) and a(€) are necessarily different, it may be pos- 
sible to use the abundance of giant arcs as a function 
of length, width, and length to width ratio simultane- 
ously in order to remove source radii dependences on 
the predictions. Alternatively, it may be the case that 
giant arc abundances are better suited to provide con- 
straints on the proper ties of high redshift galaxies (e.g., 
iBezecourt et al.lll998|) . We shall not attempt to resolve 
this question here. 

The general distinction between distortion and image 
merging arcs we have introduced also has repercussions 
for semi-anal ytic calculati o ns of lensing cross sections. 
In particular, iFedeli et al.l (|2006^l argued that one may 
obtain accurate cross sections using a simple method in 
which the eigenvalue ratio of the inverse magnification 
tensor is mapped onto the source plane, and then con- 
volved with the appropriate top hat filter for the source 
under consideration. Since such an algorithm explicitly 
ignores the possibility of image merging, it is evident that 
it will underestimate lensing cross sections whenever im- 
age merging arcs provide a non-negligible contribution 
to the total lensing cross section, as shown in Figure ITDl 
Moreover, at fixed length to width ratio the underesti- 
mate will become less severe as halo mass increases, an 
e xpectation that is fu lly consistent with the data shown 
in lFedeli et al.l (I2006D. Note that for the massive clusters 
considered in IFedeli et ail l)2006j) . we do indeed expect 
excellent agreement between the semi-analytic estimates 
and the true lensing cross section. 

As a final application of our results, we consider now 
the multiplicity ratio, the number of clusters with multi- 
ple arcs over the number of cluste rs hosting a single arc . 
This ratio was first introduced bv lGladders et aiTl)200. c !f ) 
to argue that there must be a population of super lenses: 
clusters of ordinary mass that for some unknown rea- 



son, most likely projection effects, have extraordinarily 
large cross sections relative to other clusters of compa- 
rable mass. This raises the interesting question of how 
does baryonic cooling affect the multiplicity ratio for a 
fixed cosmology. 

We estimate the impact of baryonic cooling on the mul- 
tiplicity ratio as follows: the probability P(2 + |l + ) that 
a cluster with at least one arc host two or more arcs is 
given by 



P(2+|r 



dnP(2 + \riP(n\l + ) 



(15) 



where P(2 + |/k) is the probability of observing more than 
one arc in a halo where the expected number of arcs is 
/i. P([i\l + ) is the probability that a halo hosting at least 
one arc has an expectation number of arcs /i. By Bayes's 
theorem, 

P(1+| M )P(m) 



p{n\r 



(16) 



P(l+) 

and the probability that a halo host at least one arc is 



P(l+) = J d[i P(l+|/i)P( M ). 

Putting it all together we find 

/ dix (J>(2+|/i)P(l+| M )P( M ) 



p(2+|r 



Jdfi p(i+|/i)P( M ) 



(17) 



(18) 



To compute the multiplicity ratio, we assume Poisson 
statistics for P(N\fi), the probability of a halo hosting 
N arcs given its expectation value \i. For P(/u) we use 
the fact that for a halo of mass to, [l oc a on m 7 , so 
that P(/x) cx /z _a+7_1 where a is the slope of the halo 
mass function at the lensing scale. Finally, the integrals 
of the above expression are truncated at a lower limit (Iq 
where we naively expect 0.1 < fio < 1. Assuming fj,Q is 
the expected number of giant arcs due to clusters at the 
cutoff mass m a( i or m coo i , we obtain 



(Mo) 



cool 



(W))r 



' cool (jTlcool) 

<Jad(rn acl ) 



A// 2 



1.5 



(19) 



where to obtain a numeric value we used 7 ~ 1.3 as 
appropriate for an SIS. Thus, even though cooling low- 
ers the lensing mass cutoff, it may slightly increase the 
effective cutoff in /i, and thus slightly enhance the mul- 
tiplicity ratio. Perhaps more important, however, is to 
recognize the fact that our above argument shows that 
the natural expectation value for the multiplicity ratio is 
in the range « 0.2 — 0.6 as can be directly computed from 
equation^] and that this ratio should be fairly robust 
to details of the halo p rofile and cosmology, as observed 

bv iHo fc white! Fmm . 

8. SUMMARY AND CONCLUSIONS 

Using simple analytic halo density profiles, we have in- 
vestigated how lensing cross sections depend on various 
halo parameters. In particular, we found that the lensing 
cross section has three distinct contributions correspond- 
ing to arcs formed through image distortion only, merg- 
ing of two images, and merging of three images. More- 
over, a(£) exhibits a knee at the scale which the cross 
section goes from being distortion dominated to image 
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merging dominated, and the latter component falls ex- 
ponentially fast with increasing length to width ratio. 

We then proceeded to investigate how these various 
contributions to the lensing cross section depend on halo 
properties. In particular, we found that while the image 
merging contribution to the cross section increases with 
increasing halo cllipticity, the distortion contribution re- 
mains relatively constant, except for the fact that it gets 
truncated at a smaller transition scale £ t - In other words, 
It decreases with increasing cllipticity. Wc also investi- 
gated the impact of the slope of the profile on the lensing 
cross section, and found that at fixed Einstein radius, 
steeper profiles result in more effective lenses, and more 
circular distribution of arcs about the cluster center. We 
also noted that the ratio of cusp arcs to fold arcs appears 
to be very sensitive to the slope of the density profile. 

Based on these observations, we argued that the prob- 
ability for finding arcs of a given length to width ratio 
should quickly rise with mass below some mass threshold 
M m i n , and that above this mass cut the scaling should 
become much flatter. We also discussed how the distinc- 
tion between arcs formed through image merging and 
arcs formed through distortion leads to a non-trivial de- 
pendence of the lensing cross section on the assumed 
source radius. 

We also investigated the impact of baryonic cooling on 
the lensing cross section, and found that the steepening of 
the halo density profile in response to the central baryon 
condensation enhances the lensing cross section at low £ 
by about a factor of three, and lowers the mass threshold 
M m i n for halos to become effective lenses. The latter ef- 
fect helps further increase arc abundances by about 50%. 
Moreover, cooling clusters had a significantly larger prob- 
ability of forming fold arcs than adiabatic clusters, and 



Bartelmann, M., Huss, A., Colberg, J. M., Jenkins, A., & Pearce, 

F. R. 1998, A&A, 330, 1 
Bartelmann, M., Meneghetti, M., Perrotta, F., Baccigalupi. C, & 

Moscardini, L. 2003, A&A, 409, 449 
Bartelmann, M., Steinmetz, M., & Weiss, A. 1995, A&A, 297, 1 
Bartelmann, M. & Weiss, A. 1994, A&A, 287, 1 
Bezecourt, J., Pello, R., & Soucail, G. 1998, A&A, 330, 399 
Blumcnthal, G. R., Faber, S. M., Flores, R., & Primack, J. R. 1986, 

ApJ, 301, 27 
Cooray, A. R. 1999a, ApJ, 524, 504 
— . 1999b, A&A, 341, 653 

Dalai, N., Holder, G., & Hennawi, J. F. 2004, ApJ, 609, 50 
Fedeli, C, Meneghetti, M., Bartelmann, M., Dolag, K., & 

Moscardini, L. 2006, A&A, 447, 419 
Ferland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. W., 

Kingdon, J. B., & Verner, E. M. 1998, PASP, 110, 761 
Flores, R. A., Mailer, A. H., & Primack, J. R. 2000, ApJ, 535, 555 
Gladders, M. D., Hoekstra, H., Yee, H. K. C, Hall, P. B., & 

Barrientos, L. F. 2003, ApJ, 593, 48 
Gncdin, O. Y., Kravtsov, A. V., Klypin, A. A., & Nagai, D. 2004, 

ApJ, 616, 16 

Grossman, S. A. & Narayan, R. 1988, ApJ, 324, L37 

Haardt, F. & Madau, P. 1996, ApJ, 461, 20 

Hamana, T. & Futamase, T. 1997, MNRAS, 286, L7 

Hattori, M., Watanabe, K., & Yamashita, K. 1997, A&A, 319, 764 

Hennawi, J. F., Dalai, N., Bode, P., & Ostriker, J. P. 2005, ApJ 

Submitted, astro-ph/0506171. 
Ho, S. & White, M. 2005, Astroparticle Physics, 24, 257 
Horesh, A., Ofek, E. O., Maoz, D., Bartelmann, M., Meneghetti, 

M., & Rix, H.-W. 2005, ApJ, 633, 768 
Kaufmann, R. & Straumann, N. 2000, Annalen der Physik, 9, 384 
Kazantzidis, S., Kravtsov, A. V., Zentner, A. R., Allgood, B., 

Nagai, D., & Moore, B. 2004, ApJ, 611, L73 
Keeton, C. R. 2001, ApJ, 562, 160 
Kcnnicutt, R. C. 1998, ApJ, 498, 541 

Klypin, A., Kravtsov, A. V., Bullock, J. S., & Primack, J. R. 2001, 

ApJ, 554, 903 
Kochanek, C. S. & White, M. 2001, ApJ, 559, 531 



their lensing cross section was much less dependent upon 
the line of sight projection axis than adiabatic clusters. 
Both of these effects are also explained by the steepening 
of the halo profile brought about by the contraction of the 
halo in response to the condensation of cooling baryons 
at the center. Finally, we also argued that substructures 
did not appear to significantly enhance lensing cross sec- 
tions, and that the multiplicity ratio, i.e. the number 
of lenses that exhibit multiple arcs, is only moderately 
affected by the details of the mass distribution of the 
lenses. 

Given the sensitivity of arc abundances to baryonic 
processes, it is difficult to see this observable becoming 
a tool for precision cosmology in the near future. Con- 
versely, as the range of cosmological parameters get nar- 
rower from other observations, this same sensitivity can 
be used to probe cluster properties and/or the properties 
of high redshift galaxies. Concerning this last possibil- 
ity, it seems clear that a systematic analysis similar to 
the one carried here aimed at understanding how various 
source properties affect lensing cross section would be a 
fruitful endeavor which we intend to carry out. 
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Fig. Al. — Effects of ring images on the lensing cross section o~{l). The triangles and diamonds points show a{l) as a function of L/W 
with and without ring images for an SIS profile with velocity dispersion o~ v = 10 3 km s and lens and source redshifts z; = 0.3 and 
z s = 1.5. This shows that failure to treat rings as such can result in biased lensing cross sections of circularly symmetric profiles. The solid 
line fit is our analytic fit composed of a power-law distortion component {dotted line) plus a constant contribution corresponding to the 
cross section for forming rings {dashed line). Note that the total cross section with rings consists of the cross section without rings plus a 
constant contribution corresponding to the cross section for forming rings, nR 2 . 
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APPENDIX 
LENSING ALGORITHM 

In this section, wc describe the details of our lensing algorithm, including various procedures to secure a 10% accuracy 
in the cross section and several techniques to improve the speed of the algorithm. 

Source Placement, and Image Identification and Classification 

Having defined the lens and source plane grids, we define the placement grid resolution r p via r p = R/N with N = 5 
as our default value, which is fine enough to adequately resolve the source plane region where sources touch the caustic 
of the lensing cluster. Now, for our most massive systems, the placement grid may easily have > 2000 pixels per side 
with < 1% of these pixels producing giant arcs. We therefore seek an efficient algorithm for source placement and 
image processing. 

We proceed as follows. We begin by first placing sources on all grid pixels that map onto a region where Ar/At > 
0.9£ m i n and iteratively place sources along the edges of the above region until the placement of additional source do not 
result in new arcs. We record the position, length, and width of every image generated, as well as the corresponding 
source position. 

The images of an individual sources are processed as follows. First, a list of all lit (not necessarily contiguous) image 
pixels is made. A random pixel is selected, and we move along the x and y axis until we hit the edge of the image, 
collecting all lit pixels as we go along. For each pixel collected, we then move in the orthogonal direction, and iterate 
until no new pixels are collected. We also keep track of edge pixels. The output of this procedure is thus a list of each 
of the disjoint images along with the edges of each of these images. We note that the procedure described above is 
much more efficient than the simple neighboring pixel search which is typically employed. 

Consider now a single arc. To compute the length and width of the arc we first identify the image pixel containing 
the center of the source (call it C), the pixel farthest from it (call it F\), and the pixel F^ farthest from F\. The 
length of the image is defined to be that of the circular arc defined by the points F\ C F^ plus one pixel unit, the latter 
being a geometrical correction introduced bv iPuchwein et all 1(20051 ) . Finally, we define the width of the image via 
irLW = AA where A is the area of the image, L is its length, and W its width. Our definition is motivated by the fact 
that for an ellipse, ttLW — AA. 



14 



Rozo et al. 



Circular Profiles and Formation of Rings 

There is one special case that our algorithm does not treat properly: that of sources that produce full Einstein rings, 
since for the latter A sa 2irLW where W is the width of the ring. The obvious way to correct for this is to classify 
images as arcs or rings, but this can be laborious. Instead, we have opted for a statistical approach. The basic idea 
is simple: assuming rings are formed, the width distribution of the longest arcs in the sample will be bimodal with a 
large gap of about a factor of two corresponding to the transition between disjoint arcs and rings. We can thus identify 
rings by searching for this gap, and then identifying all of the wide "arcs" as rings. Operationally, we identify all arcs 
with lengths and widths larger than some length and width cuts L cu t and W cu t as rings, where we set L cu t to be 80% 
of the length of the widest arc and a width cut W cu t to be the median width of all arcs with L < L cut . Figure [XTI 
shows the effects of rings on the lensing cross section a{£). It shows that failing to treat rings properly resulting in 
an underestimate of the giant arc cross section for high £ values since rings are being misidentificd as relatively small 
(low £) arcs. To correct for this, in this work we simply assign to all rings a length to width ratio £ =max(f) where 
max(f) is the maximum length to width ratio of all other arcs. 

Selection Function: Choosing Tangential Arcs 

Arcs can form in one of two ways: through strong distortion along the tangential direction, or through the merging 
of different images of a single source. This is true for both tangential and radial arcs. Therefore, the cross section 
a(£) generally consists of four distinct contributions. However, we have chosen to focus on tangential arcs only, 
since tangential arcs are much more prominent and easily identifiable than radial arcs. Observationally, this has the 
advantage that the purity and completeness of a tangential arc sample will be higher than that of a radial arc sample. 
It is therefore important to identify and remove radial arcs when computing the cross section for tangential arcs. 

We have chosen to identify radial arcs through a statistical procedure similar to the one used to identify rings. 
In particular, we exploit the fact that radial arcs tend to be not just magnified along the radial direction, but also 
demagnified along the tangential direction. In other words, the radial arcs create a distinct population of images 
branching out from the main population and reaching the low L and low W part of the plane. We therefore tag all 
arcs with L < L cut and width W < W CM t as radial arcs. The cuts are defined through the following algorithm: first we 
find the median width of all arcs above the minimum length to width cut £ m in, which is our first estimate for W cu t- 
We then lower W cu t in steps of 0.5r p where r p is the placement grid resolution until the increase in the number of arcs 
with W > W cu t is less than 5%. The length cut L cut is defined as L cut = W cu t£ m i n where £ m in is the minimum length 
to width ratio of arcs considered. We found these cuts cleanly separate the tangential and radial branches. 

VALIDITY OF THE SMALL SOURCE APPROXIMATION 

One of the key steps in our lensing finding algorithm is selecting all source placement grid pixels where the eigenvalue 
ratio Ar/At is larger than 90% of the minimum length to width ratio considered. This selection is based on the 
expectation that for infinitely small sources, the eigenvalue ratio Ar/At will be precisely the length to width ratio £ 
of the resulting image. A natural question that arises then is to what extent is this small source approximation valid. 

Figurc lBTI illustrates the relation between L/W and Ar/At for an SIE model with a v = 10 3 km s _1 and q = 0.8. Each 
point in the figure represents a source that does not touch the lens's caustic, while sources that touch the tangential 
caustic are shown as diamonds. This shows that the agreement between Ar/At and £ is excellent provided the source 
in question does not touch the lens's caustic. This str ongly suggests that the generalization of the eigenvalue ratio 
Ar/At estimate for £ to ellip tical sources proposed by iKeetonl 1)20011) should be accurate for distortion arcs, as was 
indeed found bv lFedeli et all <|2006|) . 

FITTING FUNCTION FOR LENSING CROSS SECTIONS 

In this section, we develop an analytic fitting function to the lensing cross sections useful for both analytic and 
numerical models of galaxy clusters. The idea is based on the realization that at low £, the lensing cross section curve 
is dominated by distortion arcs is well fit by a power law, whereas at large £ ratios, the cross section is dominated by 
image merging arcs, and is better fit by an exponential fall off. So, the total cross section a{£) can be expressed as a 
sum of these two components: 

— O~distortion(£*) ~t~ Emerging (^) (^-0 



where 

0-distortion(£) = \ °"° ~ (^^) 




if £ >£ 



d^raax > 



and 



a c eiq>[-(£-£ t )/£ s ] if £ > £ t , 



In the above expressions, £dmax is t ne maximum length to width ratio of distortion arcs, and £ t is the minimum length 
to width ratio of image merging arcs. The subscript c denotes 'caustic', in that image merging sources are those that 
touch the lens's caustic. 
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Fig. Bl. — I vs. eigenvalue ratio. The plot is constructed as follows: we begin with an SIE profile with a velocity dispersion a v = 
10 3 km s _1 and an axis ratio q = 0.8. For each source we place on the source plane, we find its largest arc, and the eigenvalue ratio Ar/At 
at the center of said arc. To keep the x-axis small, we set Ar/At = 90 for sources in which Ar/At > 90. Dots represent sources which 
do not touch the tangential caustic, whereas diamonds correspond to sources that do touch the tangential caustic. It is evident from the 
figure that the linear approximation £ = A^/Ax, shown here with the solid line, is a good one even for strongly distorted images, but that 
the agreement fails if the source touches a caustic. This reflects the fact that the length of these arcs is set by the merging of two images. 

While £d,max tends to be difficult to resolve 13 , we may simply set £d,max = 4 in the above equation. This results in 
the fitting function: 



{a c exp[-{£-£ t )/£ s } ii£>£ t 

Our ansatz in eauation lC4l providcs an excellent fit to the data whenever the distortion and image merging contributions 
can be cleanly separated. 14 All fits were done with logarithmic sampling on £, and assuming an arbitrary 10% error 
bar at each point. 15 Furthermore, the x 2 surface typically exhibits only one deep minima, implying that degeneracies 
between the various parameters do not exist when the two components of the lensing cross section are clearly distinct. 

Note that every parameter has a concrete physical meaning, but these parameters are clearly correlated. For example, 
we expect the parameters £t and £ s to be correlated, since £ s is a measure of how quickly the length to width ratio of 
an arc decreases as the source moves away from the lens's caustic, and £t is the minimum length to width ratio such 
a source can have while still touching the caustic. Another interesting correlation is between the parameter 4 and 4 
where 4 is defined via 

t c = <to{£o/£*) v - (C5) 

Thus, 4 is the length to width ratio at whic h th e image merging cross contribution to the cross section is as large 
as the image distortion contribution. Figure IC1I shows that a clear and tight correlation exists between £ t and £ s 
(dashed line) and between £ t and 4 (solid line) for both analytical and simulated cluster data spanning a large range 
of ellipticities, halo profiles, and Einstein radii. The best fit lines to the 4 _ 4 and 4 _ 4 correlations are given by 

4=0.748 + 0.4094, (C6) 
4 = 0.260 + 0.6204- (C7) 

Note that in Figure IUT1 the 4 values have been displaced upwards by five units for illustration purposes. 

The presence of strong correlations among fit parameters suggests that five parameters are not needed to adequately 
fit the data; three parameters will suffice, namely a®, v, and 4- It turns out that there is a tight correlation between v 

13 We often found that id, max &t> and max tends to be difficult to resolve since the cross section at these scales is dominated by 
the image merging contribution. 

14 This is essentially all cases we considered with the exception of most (3 = 0.6 models, and numerical clusters where the presence of 
substructures severely distort the shape of the lensing cross section. 

15 These do not represent real errors, which will in reality be strongly correlated, and are simply used to define a best fit model. 
Consequently, the \ 2 f° r these fits cannot be used to judge whether the fits are statistically acceptable or not. We choose our particular 
fitting scheme since in the absence of a clear way to assign error bars to our measurements, simplicity was the next obvious criteria for 
choosing a fitting algorithm. 
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Fig. CI. — Correlations between the fitting parameters of equation IC4I In addition to the correlation between it and £ s shown by a 
dashed line, there is an additional correlation between It and the parameter £* defined through oo(lo / = cr c . This latter correlation is 
shown with the solid line. We have used a least square fitting to define the correlation parameters, and have considered only data in which 
the cross section function exhibited a clear transition scale, with the added provision that the image merging contribution to the cross 
section by dominated by fold arcs. Note that the £* values have been displaced upwards by five units for illustration purposes. Also shown 
above are the best fit values obtained from x, y, and z projections of simulated clusters including cooling and star formations, divided into 
clusters with a relatively isolated central component (triangles) and clusters with large amounts of substructure near the central density 
peak (crosses). Note that the best fit relation obtained from our analytic profiles appears to be satisfied by numerical clusters as well. 

and i t as well, further suggesting that lensing cross section curves o~(£) may, in general, be characterized through only 
two parameters, an amplitude <7o and a transition scale if This is intriguing since the minimal number of parameters 
to describe a halo are three: an Einstein radius, an ellipticity, and slope profile /3, which raises the interesting question 
of whether it is possible for two different halos to have the exact same lensing cross section function cr(t). 



